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ABCTRAGT ^ 



The O , ^ tracker system of optimal smoothing was investi- 



gated in a track-while-scan radar by using three dimensional simulation of 
various interceptor-target profiles. Random noise was put on all parameters 
which were not positively defined. This included the target's motion. 



The computer program used is designed for carrying on the analysis 
with other systems so that their performance can be compared with the 



which should provide a suitable basis for evaluation of any chosen 
alternate system. 

The parameters used in the analysis are similar to those which 
will be used in the Hughes Aircraft Co. developed AWG-9, Phoenix, 
Aircraft Missile Control System designed for the F-lllB aircraft. 




tracker. A number of representative profiles are listed 
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Introduction 



The ^ ^ sampled data system has been the most often 

proposed means of performing the task of optimum smoothing in a 
track-while-scan radar. This method is particularly attractive since 
it is very simple. However, in the light of the high performance intercept 
problems which may be encountered today, especially by a high perfor- 
mance aircraft, it was felt advisable to investigate the capabilities of 
such a smoothing system in the environm^-nt of air-to-air intercepts. 

In this analysis the third order system or 
tracker was considered by using a three dimensional simulation of various 
interceptor-target profiles. Random noise was placed on all the para- 
meters which were not positively defined. This included the motion of 
the target. 

The computer program used was designed for carrying on the analysis 
with other systems so that their performance could be compared with that 
of the 0< , tracker. 

The parameters used in the analysis were similar to those which will 
be used in the Hughes Aircraft Company developed AWG-9, Phoenix, 
Aircraft Missile Control System. This system is currently being developed 
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for use in the F-lllB aircraft. 
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2. Basic System Description 

The system considered in the analysis was a pulse doppler radar 
operating in a track-while-scan mode, TWS . That is, the radar was 
able to obtain simultaneously at discrete sampling intervals, called 
frame times, the target's range, range rate, and bearing, in azimuth and 
elevation, relative to an interceptor’s antenna. At the same time the 
system was also maintaining a track file on a number of other targets 
at locations spaced over its scan volume . 

The heart of a system such as this is the so-called "track-while- 
scan loop. " Fig. 1 is a basic block diagram of this loop. 

Prior to entering the TWS loop the raw A to D converted parameters 
are taken from the output of the radar, and the interceptor's inertial 
reference platform and air data computer. This information is pre- 
processed to remove any redundanciescaused by picking up the same 
target on two or more antenna sweeps, called bars. Also the north, 
east, down direction cosines are computed. These are the direction 
cosines for a north, east, down axis coordinate system which is fixed 
to the interceptor. The axes are defined by having the north axis pointing 
toward true north. The east axis results from the cross product of the 
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down and north axis. The down axis points towards the earth's local 
verUcal, and is positive downward. Thus an interceptor flying north 
would have the north axis along the line of sight, the east axis out the 
right wing and the down axis vertically downward . 

The data now enters the "track-while-ican loop." It is first 
correlated by having the individual observations checked against previous 
tracks on file in memory. This may result in several observations being 
linked with one track. These are then associated by having the observations 
paired by a logic routine with the track which it most closely matches . The 
next stage initiates new tracks or deletes old tracks , which no longer 
associate after a quality counter, which rates each track, degenerates to 
a point such that the track no longer warrants being maintained on file. 

The next section and the one which is investigated in this paper is 
that of smoothing and prediction. The ability of the system to perform all 
previous and following tasks rests on how well the track is smoothed, i.e. , 
the effect of noise removed, and how well it predicts the target’s future 
motion. This section must be designed to respond quickly to target 
maneuvers while filtering out the false excursions of the target caused 
by noise so that future predictions will be as accurate as possible. 
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The gating section is next in the loop. It is this section which 
determines the possible limits of target motion between contacts. It allows 
for target maneuvers between scans plus approximately 3cr noise. It further 
provides for expansion in the event of a missed contact. A maximum gate 
size is imposed which corresponds to a volume which would contain any 
possible target maneuver that could occur up to the time the track would 
normally be eliminated by the quality counter due to missed contacts . 

The individual track parameters are stored in the track file section 
to keep it updated. These along with the track locations are used to 
determine target priorities . The track parameters are also used to determine 
which observations correlate with the tracks during the correlation phase 
at the beginning of the loop. The above computations are performed on all 
tracks at the end of each frame time . 

Additional functions which result from the TWS loop computations 
are the antenna scan center calculations to keep the pattern centroid 
located in the center of the incoming raids , the steering information to the 
pilot to allow him to maneuver as necessary to keep the radar within its 
gimbal limits, prelaunch information to the missile so that its seeker head 
is oriented in the proper position at the beginning of its semi-active phase. 
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and midcourse maneuver corrections to the missile during flight. 

It is readily apparent that the system must be optimized to the highest 
possible degree in light of the high relative speeds encountered and the 
number of targets to be tracked. At the same time it is extremely important 
that computations remain relatively simple to maintain reliability and to 
keep computational time and storage space to a minimum. 

The ability of the system to meet these criteria lies primarily in the 
smoothing, prediction and gating sections. It is the former two which are 
considered in this paper. 
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3. System Characteristics 



The values of the system performance parameters were derived 
primarily from the Hughes SYSTEM EVALUATION ^ND ANALYSIS REPORT , 
SEAR, /l^. Where possible these parameters have been crossed checked 
with other company internal documents , or by conversations directly with 
company engineers. The following parameters and assumptions are those 
which were used in the simulation. 

It was assumed that the interceptor would have an on station Mach 
number of 0.7 at above 36,000 feet. The target's speed was believed to 
be between Mach 0.4 and 3.0 and its altitude could be between sea level 
and 80 ,000 feet . 

The general characteristics of the radar were those of a high PRF 
pulse doppler using linear frequency modulation. 

The detection range was 165 n.m. This was defined as the range 
where S/N was one using the classic radar range equation. It was based 
on a five square meter target. The radar's probability of detection curve 
is shown in Fig . 2 . 

The effect of clutter on detection probabilities was not considered 
since a large amount of clutter rejection is an inherent advantage of a high 
PRF pulse doppler radar. This is particularly true over water. In addition 
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the system will also have clutter rejection circuitry. Due to the high PRF 
certain relative velocities may cause periodic loss of signal due to 
"eclipising" , i.e. , target echo returning during pulse transmission. 

Multiple PRF rates will normally remove this possibility. Thus eclipising 
was not considered to effect detection probability. 

In the radar, range is measured by taking the difference between 
two filters in a doppler filter bank that are "rung" on successive FM phases. 
The spacing of the filters corresponds to five n.m. Since there are 30 
filters the range can be measured to 150 n.m. The signal was considered 
to originate uniformly in plus or minus one-half a filter and terminate 
normally distributed with a standard deviation of one-third a filter. 

This is shown in Fig. 3. 
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Hughes stated that a joint distribution describing these two distributions 
yielded a normal distribution of zero mean and standard deviation of 
1.37 n.m. It was felt that this value was too optimistic and indicated that 
the distributions were quite dependent , i.e. , the covariance was large. 

The value used for the range standard deviation in this analysis was 
2.5 n.m. 

Range rate is measured by frequency shift as determined by a doppler 
filter bank. The filters are spaced 150 cps apart. Thus at X band using 
the equation given by Skolnik /_ 1 _/: 



The total number of filters is 768 giving the ability to measure range 
rates of up to approximately 3250 kts. 

It was determined by Hughes that the distribution of the range rate 
error was normal with zero mean and a standard deviation of two filters , 
8.72 kts . 



V = ^d ^ 
r 1 no 



103 



Where: f^ = doppler frequency shift. 




wavelength in cm . , 3 cm at X band . 



V = range rate in knots . 
r 
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The scan pattern of the antenna in the TWS mode is either two bar 
+ 40 deg. or four bar + 20 deg. The beam separation is 1.32 deg. and the 
scan rate is 40 deg/sec. It was determined from various sources that the 
average antenna azimuth error could be considered normally distributed with 
zero mean and standard deviation of 0.745 deg. Similarly the elevation 
error had zero mean and standard deviation of 0.86 deg. These values 
include radome , potentiometer output, and pick off errors. The additional 
error in elevation is due primarily to the introduction of redundancies caused 
by picking up the same target on two separate scans . 

The smoothing and prediction system which has been proposed by 



Hughes is the 




, 'X' tracker type. The characteristics of this 



system are considered in the following section. 
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4. 



Smoothing System Discussion 



Since computational time and storage space available for smoothing 
and prediction is quite limited some of the more sophisticated smoothing 
systems would not at first appear to be suitable for use in the system 
described. The requirements of the system are reduced somewhat, however, 
by the fact that the target’s capabilities are relatively predictable, but of 
course, entirely random. For these reasons the most attractive choice is 



be acceptable would reduce to a scheme which would use a similar smoothing 
constant gain technique in any case. 



tracker to allow for the smoothing of acceleration. As given by Benedict 
and Bordner / 3/ the variance reduction ratio is optimized by having the 



an 




tracker type of system. Almost any other system that would 



The 




, tracker system is an extension of the , 



velocity constant, , defined as : 



z 




Similarly it was shown by Simpson l_\J that when acceleration was 
considered the optimum value of the acceleration smoothing constant, Y” , 



became: 
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Thus if it is possible to define the smoothing equations 
become unique . 

The proposed equations used for the smoothing and prediction 



are: 



A. Range: 






'R-. 






Where: 



o if an observation 



was received. 

and v^ observation was missed. 

T-Frame Time = 2+0.1 sec. 

c^= Number of frame times since last value of range rate was 
received . 

The subscripts are: 
n - Frame number. 

P'= Predicted . 

O = Observed . 

) . • ) * ) » > < I » > t » 

'» » k 

S)nooilied,. ’ J ^ 
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The value of for range smoothing as proposed by Hughes was 
0.45. Since the standard deviation of range rate was relatively small 
no smoothing of range rate was necessary so that no ^ , or 
were required in range smoothing. 

B. Direction Cosines: 



JVn-” '■T_)V(n.O * 

(jV.„- -JVp^ 



Where: 

J\-.‘ A,., if an observation was received. 

A.- if an observation was missed. 

A = general direction cosine being smoothed . 
For the North and East direction cosines, N, E, Hughes proposed 
the smoothing constant values of: 

^ 0.43 

0.120 

0.008. 
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For the down direction cosine, D, the proposed values were: 

^ - 0.35 
=0.075 
0.0035. 

A Z- transform analysis of the o<^ , ^ Tracker was performed by 
Sklansky /sy. Basically the system is stable and for a non-maneuvering 
target there would be no steady state error for between zero and one. 
When is equal to one no smoothing would be performed and raw data 
would be the output. In a high noise non-maneuvering target condition 
the smoothing constant would be small and the system response would 
be slow. For increasingly lower noise the smoothing constant would 
approach one in the limit. For a maneuvering target the smoothing 
response must be increased so that a larger value of is required. 

The obvious disadvantage of the Tracker is that there is 

not a single optimum value of the smoothing constant for all the situations 
that can be encountered. 

The advantage of the system again lies in its simplicity and the 
small number of computations and amount of storage space required. 

Also, since and ^ are defined uniquely when is determined the 
problem of fixing the constants is reduced somewhat. 
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5 . Simulation 

The CDC 1604 computer was used to implement the three dimensional 
simulation of a single target-interceptor problem. This program was named 
Three D. The primary aims of the simulation were to provide realistic, 
flexible and accurate target-interceptor geometries . The pseudo random 
noise generator, sub routine RNDEV, was used to provide gaussian noise 
on all parameters that were not positively defined . 

A list of symbols and a description of the use of the programs used 
are contained in Appendices A and B. 

The resulting program allowed the target to be positioned initially 
at any location in the northern half of a north, east, down coordinate 
system which was ground stabilized at the interceptor's initial location 
at problem time zero. The target's initial coordinates, maneuver para- 
meters, and maneuver times were then read in by data cards along with 
the interceptor's heading, altitude, and Mach number. The interceptor's 
velocity was computed in the program by using the following relation: 

Vfpg = Sonic Velocity fpg x Mach Number. 

The variation of the speed of sound with altitude is shown in Fig. 4. 
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Fig. 4, 

Speed of Sound Versus Altitude 




No provision was included to maneuver the interceptor since one of 
the basic assumptions of the system was that the interceptor was continually 
tracking several targets in widely spaced locations in the radar s scan 
pattern. It was thus not possible for the pilot to place an individual target 
in the center of the scan. Also no random noise was used to vary the inter" 
ceptor's position since the same relative effect was gained by randomly 
varying the target's velocity components. 
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The problem time was taken to be 300 seconds . This allowed a target 
at Mach 1.5 above 36,000 feet to travel slightly over 70 n.m. in the course 
of the problem. Since the probability of acquiring and tracking a target 
beyond 100 n.m. was very low, 100 miles could be considered the range 
limit in setting up the problem profile . 

The target's velocity was also determined by the previously mentioned 
relations. No provision was made to accelerate the target, but the target 
would decelerate at the start of its dive. Its final speed was a function 
of altitude below 10,000 feet. The minimum speed was Mach 1.2, 1340 fps. 
at sea level. The relations used are shown in Fig. 5. The target was 
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Fig. 5. 

Final Target Velocity Solow 10,0C0 Fe^at ‘ 
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programmed to make a 180 degree turn at a programmed time and horizontal 

turn g load. It would also make a 60 degree dive and pull out at a pro- 
grammed time, final altitude, positive and negative g factor. 

For convenience, distances were read in, in nautical miles, velocities 
in Mach number, and headings in degrees. The program then converted the 
values to compatible units of ft. , ft. /sec. , and radians. 

The motion of the target was computed by piece-wise-linear integration. 
Random accelerations having a standard deviation of one-half g for two second 
periods were imposed. The target's motion was broken down into velocity, 

VT, along its path of flight, heading angle, AT, and pitch angle, BT. Each 
component was then related to the distance traveled in each of the coordinate 
axes in the iteration time and combined with the previous position to give 
the new target location , XT , YT , ZT . 

Every two seconds of problem time the target's coordinates were taken, 
compared with the interceptor's coordinates XO, YO, ZO and the actual 
antenna azimuth, ATA, elevation, EPSI, range, R, and range rate, RDOT 
were computed. If the target had passed out of the limits of the radar beam,+ 

65 degs . in azimuth and elevation the problem stopped. It should be noted 
that the antenna does not have a scan volume this large in the track while 
scan mode , but it does have the capability of positioning its scan anywhere 
inside these limits,. .. _ .... .... .. .. 
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The random noise generator was then used to put noise on the radar 
parameters, ATAE, EPSIE, RE, and RDOTE, and the apparent target location 
having coordinates XTE, YTE, and ZTE, was determined. Since a relatively 
small number of samples were used the mean, variance and independence 
of the noise for each of the parameters were checked using an autocorrelation 
routine. The results of analysis are shown in Table I. 

The actual and noisy direction cosines for the three axfes' were then 
computed. These were, north direction cosine, DIRN , and DIRNE, east 
direction cosine, DIRE and DIREE, and down direction cosine, BIRD, and 
DIRDE . 

Knowing the probability of detection of the radar, RNDEV was used 
to determine a zero-one probability of detection for each frame time . This 
was done by cross plotting the absolute value of the deviation required, 
versus range, corresponding to the probability of detection curve. Fig. 2. 

This was then made piece wise linear and implemented by using IF statements, 
Thus for each actual range a maximum value of deviation for detection could 
be determined, if the absolute value of DEV obtained from RNDEV exceeded 
this value no return for that frame time was assumed. Fig. 6 illustrates the 
method used. 
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TABLE I 



AUTOCORRELATION ANALYSIS RESULTS OF PROGRAM 

THREE D 
(151 SAMPLES) 





Mean 


Variance 


Largest Value of 
Autocorrelation 


RN 1 
Target 

Acceleration 

Noise 


0.0112 


.772 


0.1787 


RN 2 

Target 

Heading 

Acceleration 

Noise 


-0.001 


1,0636 


0.1706 


RN 3 

Target 

Pitch 

Acceleration 

Noise 


-0.1032 


0.9585 


0.2358 


ATA 

Antenna 

Azimuth 

Noise 


0.0785 


0.9612 


0.1717 


EPSI 

Antenna 

Elevation 

Noise 


0.6308 


1.0433 


0.1576 


R 

Radar 

Range 

Noise 


0.1386 


0.9938 


0.1721 


RDOT 

Radar 

Range 

Rate 

Noise . . . .. 


0.1114 

♦ * 1 I 


0.9780 


-0.1528 


PD : ■ : 

Radar ’ ’ ‘ ’ * 

Probability 
of Detection 
Noise 


• l* 

• ) 

* * • • • > > 

-0.0564 


• • > 

) « * 

0.9914 


■ * 

-0.1468 
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In the actual system this would be the status of the individual track 
in the TWS loop after being correlated and associated. 

A quality counter, NQ, was devised for the program which incremented 
in steps of two for each received return to a maximum of ten. For each 
missed return it degenerated by one. When the quality counter reached 
zero the track was dropped and no smoothing was performed until the next 
received return. 

The particulars of the smoothing and prediction section were as 
described earlier. 

The primary quality indicator for the system was the three dimensional 
r ms smoothed difference from the actual position. Since high relative 
velocities between the missile and the target are considered very likely 
in the active phase of the missile flight it is necessary that the missile be 
positioned accurately at the beginning of this period so that its maneuvers 
may be kept to a minimum. In order to assist in achieving this a desired 
rms smoothing error limitation of approximately one-half mile inside 50 miles 
was imposed. On the output curves points were plotted indicating these 
limits. For range constant lines of points at plus or minus one-half mile 
were plotted. For azimuth and elevation points were plotted as a function 
of range using the relation: 
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If all the parameters were smoothed to the desired value the three dimensional 
rms error, CIRSMO , could be as high as 0.875n.m. Thus on the curves of 
CIRSMO a line of points having this value were plotted. 

The program print out contained the run number, problem time, actual 
range, rms three dimensional actual error, CIRERR, rms three dimensional 
smoothed error , CIRSMO , a one or zero for detection or no detection , ND , 
and the status of the quality counter, NQ. When the quality counter was 
zero no contact was assumed and all nines were printed out. 

To observe the results of the smoothing method the sub-routine 
CALL DRAW was used to plot the interceptor and target positions in the 
north-east, N-E, plane, looking down from above, and the north-down, 

N-D, plane, looking in from the east axis. The sub^OUtlhjp was also used 
to plot the difference between the actual and smoothed values of range 
and the direction cosines, which were plotted as angular differences, and 
the values of the three dimensional rms smoothed error. 
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6. Results of the Analysis of the Proposed System 

Table II contains the important details of the problem profiles used 
in the analysis. The results of run lA are shown in Table III and Figures 
7 to 15. The results of runs 2A and 3A are shown in Appendix D. 

Three basic profiles |were considered. These were identified by the 
initial position and heading of the target with respect to the interceptor. 

For run 1 the target was positioned at 100 n.m. and was closing the 
interceptor. For run 2 the target was positioned at 10 n.m. and was 
moving away from the interceptor. In run 3 the target was positioned at 
approximately 60 n.m. at a 30 degree angle from the interceptor and was 
headed on a cross track with respect to it. The actual profiles flown by 
the target and the interceptor could be varied in a myriad of ways , however 
in reaching the conclusions that follow only those flight paths shown in 
Table II were used. 

The results obtained for the system, as proposed by Hughes, indicated 
that the smoothed tracks were still quite noisy. The direction cosine smoothing 
results obtained by analyzing the figures of these values compared to the 
desired limits indicated that the direction cosines were being smoothed inside 
the limits rather well. Some degradation was noted in the cross track cases 

• » <*•♦ '» > • » I * I 'V • ** 

) I V ♦ » * <•) I > •) 

> • > * » <( 

•• I l> ' t J W .> ,* V • * 



'25'bf lOQ 



')f* no ■>*•> 




i c 

t f 

( f 



•4« •t# 

• • • 

• • • • • 

\ * ‘ 

• •• ••• ««• 



w 



% 



H 



CO 

t—t 

CO 

< 

< 



Q 

w 

CO 

!=> 

CO 

w 

I— I 

O 

cc 

cu 

H 

o 

3 

Pi-. 




» .too 



0 ' O *.> 'too 

■> 1 

J TO > 0 

T -) 

’ T t T ■) ■) 




» J c 

• < 

• 5 . 



26 of 100 



f f ( f f I < 

f f f 

r r f f t 
r r < 

r f I f f i i 



f f t < f t t 
( < 

f ( r ( 

f r 

1 f ( f 



( f f iff 

< f 

r f < 

f r 

» f < r f f 



{ < i f f 

r f f 

f f f < 

f r f < i 

(ft* r 



f < f f f 

f r f 

t f f f 

r r * 

f f i i f f i f 



r f f 

r « f 

f f c 

f « f 

f f f 



* * ¥ 

( 

« » 



• ft • 



4 1V «* 



■ » 4* 4* 



4* 



SMOOTHING RESULTS 



RUN lA 

T = TIME, R NM = ACTUAL RANGE, C I RERR = 3d I MENS I ONAL RMS ACTUAL ERROR, 
CIRSMO-3 DiMhNS'iONAL PMS SMOOTHED !. ERROR , ND-t DETECTION, NDsO NO 
DETECTION, *N0 = QUALITY'CCUNTRR, WHRN "NQ = 0 TARGET CONSIDERED LOST AND 
ALL NINES ARE f^f^'lNTEo; * 



T 


R NM 


CIRfiRR 


CIRSMO 


ND 


NO 


f 


100.5000 


999.9999 


999.9999 


99 


0 


2. 


99.4871 


999,9999 


999.9999 


99 


0 


A. 


98.4809 


999.9999 


999.9999 


99 


0 


6. 


97.4805 


2.8089 


2.8089 


1 


2 


8, 


96,4819 


.3501 


2.7924 


0 


1 


10, 


95,4829 


999,9999 


999,9999 


99 


0 


12. 


94 ,48o8 


999,9999 


999,9999 


99 


0 


14. 


93.4664 


999,9999 


999,9999 


99 
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16, 


92,4470 


999.9999 


999.9999 


99 


0 


18, 


91,4370 


999.9999 


999.9999 


99 
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20. 


90 .4195 


999.9999 


999.9999 


99 
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22, 


89.3964 


999.9999 


999.9999 


99 
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24. 


88.3717 


999.9999 


999.9999 


99 
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26. 


87,35q7 


999.9999 


999.9999 


99 
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28. 


86,3438 


2.1180 


2.1180 
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2 


3o , 


85,3354 


1.6702 


2.0969 


0 
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32, 


84,3289 


999.9999 


999.9999 


99 


0 


34, 


83,3276 


3.94q3 


3, 94q3 


1 
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36. 


82 ,32q9 


2,5368 


3.9162 


0 
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38, 


81.3087 


[7825 


2.4350 


1 


3 


4o . 


80,2970 


3,2?25 


2.32q4 


0 
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42, 


79,2848 


.9345 


2.2242 


0 
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44, 


78,2667 


999.9999 


999.9999 


99 


0 


46, 


77.2543 


2.7971 


2.7971 


1 
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48, 


76.2567 


2.4162 


2.8147 


0 
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50, 


75,2688 


999.9999 


999.9999 


99 
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52, 


74.2774 


999,9999 


999.9999 


99 
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54. 


73,2785 


999.9999 


999.9999 


99 
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56. 


72,2763 


.7391 


.7391 
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58. 


71,2801 


3.59o5 


1.1076 


1 
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6o , 


70.2918 


2.9569 


1.2629 


0 
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62, 


69,2966 


1.8983 


1.4499 


0 
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64, 


68.2990 


2.9966 


1.5692 


1 
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66. 


67.3133 


,3755 


.7502 


1 


6 


68, 


66,3347 


3.4902 


.76q4 


0 
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70. 


65.3566 


2.4328 


.7920 


0 
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72. 


64,3877 


2.9470 


1.6471 


1 


6 


74. 


63.4255 


4 ,42l8 


1.7263 


0 
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62.4593 


4 , 4oo5 


1.8257 


0 
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78. 


61,4952 


.8885 


1.9355 


0 
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80 . 


60,5331 


.8715 


• 1.1710 
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82. 


59.5727 


3.2183 


1.201'2 
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58,6241 


6.4126 


1.2271 
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3 
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57.6929 


1.9297 


1.2447 
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2 
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94. 


53.9813 


4.2629 


1.7212 


0 


1 


96, 


53.0444 


999.9999 


999.9999 


99 


0 


98. 


52,1038 


2.9109 


2 .9109 
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100 , 


51,1811 


2.3925 


2.8871 
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5q ,2683 


1.2764 


2.1161 
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104. 


49.3610 


1.1716 


2.1280 


O 


2 


106 . 


48.4660 


1.3300 


2.1478 


0 


1 


108. 


47.5685 


2.5272 


.6311 


1 


3 


110 . 


46,6623 


1.6322 


.7560 


0 


2 


112. 


45.7441 


1.1995 


.6104 


1 
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42’,92S'2r 
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1.8194 
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l2o . 


41.9624 


1.6084 


.9878 


1 


ID 


122. 


10.9986 


3.9664 


1 , 5656 


1 


10 


124. 


4o.o442 


3.5859 


.9143 


1 


10 


126. 


39.0998 


1.04Q4 


.8934 


1 


10 


128. 


38.1587 


2,5277 


.9466 


0 


9 


130 . 


37.2168 


4.0934 


1.6227 


1 


10 


132. 


36,2832 


3,0987 


1.6066 


0 


19 


134. 


35,3557 


3.7174 


1.5809 


0 


8 


136. 


34.4268 


i.4977 


.3722 


1 


.10 


138. 


33.4964 


4.4075 


.4850 


0 


•■f9 


140 . 


32,5567 


1.4627 


.6428 


1 


10 


142. 


31.6179 


3.3400 


1.8072 


1 


10 


144, 


30 .6886 


4.5131 


1.1939 


1 


10 


146. 


29,7558 


2.1345 


1.4976 


1 


10 


148. 


28,8093 


.9501 


1.1917 


1 


10 


150 . 


27.8563 


1.1226 


,8273 


1 


10 


152. 


26.8896 


.7536 


. 6589 


1 


10 


154, 


25.9816 


3.1620 


.8572 


0 


9 


156. 


25,1572 


5.4788 


2.5853 


1 


10 


158. 


24,4380 


4.1105, 


.9185 


1 


10 


160 . 


23.8445 


1.8564 


.6228 


1 


10 


162. 


23,3933 


3.3222 


1.5404 


1 


10 


164. 


23.0952 


5.0620 


1.9082 


1 


10 


166 . 


22.9552 


.9599 


1.79o9 


1 


10 


168 . 


22.97o7 


5.6214 


2.4656 


0 


9 


170, 


.23.1315 


3l3o 


1.8425 


1 


10 


172. 


23,4222 


i.0723 


1.8749 


1 


10 


174, 


23,8221 


1.1368 


2,6427 


0 


9 


176. 


24.3fl8o 


4.4560 


.2192 


1 


10 


178. 


24.8556 


3.1212 


1.1249 


1 


10 


180 . 


25.4395 


4.2445 


2.2379 


1 


10 


182. 


26.0370 


.5878 


l.l2o6 


1 


10 


184. 


26.6251 


2.1740 


1.1881 


0 


9 


186. 


27.1849 


2.3486 


1.2085 


1 


10 


188. 


27.6986 


2.4025 


1.8340 


1 


10 


190. 


28.1555 


.6367 


2.0981 


0 


9 


192. 


28.6017 


1.5670 


1.8206 


1 


10 


194 , 


29,0530 


.7293 


.8247 


1 


10 


196. 


29.5092 


.3101 


.3816 


1 


10 


198. 


29.9700 


2.6407 


.3614 


0 


9 


200 . 


3 0 .4338 


.8438 


.4553 


1 


10 


202. 


30,9019 


2.0573 


1.0997 


1 


10 


204, 


31.3732 


1.6566 


1.2539 


1 


10 


206 . 


31.8487 


2.2964 


1.2489 


0 


9 


208. 


32.3273 


. /.77o7 


2,3333 


1 


10 


210. 


32.3C86 


:'3,1324. 


■ 2 ;. 6509 ■ 


1 


10 


2l2. 


33.2933 


: 1 , (j 363 


■ 1: 045l 


1 


1 0 


214, 


3 3 . 7 6 0 6 


i;5452 ' 


1.1828 


1 


10 


216 , 


34.2715 


1.3o69 


1.16q1 


0 


9 


2l8. 


34.7654 


.5480 


1.1220 


0 


8 


220 . 


35,26q6 


1.0041 


.3621 


1 


10 


222. 


35.7582 


2.1422 


.74x9 


1 


10 


224, 


36.2585 


1.3269 


1.00q7 


1 


10 
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TABLE III 



226, 


36 . 7614 


2.1929 


. 6306 


1 


10 


228, 


37.2662 


.4203 


.5o7o 


1 


10 


230 . 


37,7737 


1 ,4o2l 


,8871 


1 




.232, 


. .3a, 2827^ 


.2.220 9 


.6915 


0 


. 9 


;234.,: : 


38,7943j 


\.2,7^67 


: 1.6487 


4 

•L 


10 


'236; 


' •Z9:Z-(iC9 ^ 


''"'3.72 72"'’ 


" .7716 


X 


IQ 


238, 


39.8210 


.9506 


.5839 


1 


10 


240 , 


4o .3378 


.8783 


.7262 


1 


1-0 


242, 


4o , 8568 


2.7186 


1.1963 


1 


10 


244, 


■41,3762 


6 , 0 221 


1.2214 


0 


9 


246, 


■41.8975 


1.1427 


1.2499 


0 


6 


248, 


42.42o3 


3.5265 


1.2827 


0 


7 


250 , 


12.9449 


.9221 


1.0891 


1 


9 


252 , 


^3.47o4 


2,2873 


.5461 


1 


10 


254, 


43.9969 


1.0963 


.2977 


1 


10 


256. 


4^.5255 


3.24i6 


1.3572 


1 


10 


258, 


45,0546 


4,o95l 


l.l9o9 


1 


10 


26o , 


45,5847 


4.3770 


1.2842 


0 


9 


262, 


46.1162 


3.1601 


1.8915 


1 


10 


264. 


46,6486 


2.8110 


2.2411 


1 


10 


266, 


47.1817 


3,2857 


,3o59 


1 


10 


268. 


47.7164 


5.9780 


2.6170 


1 


10 


27o. 


48.2523 


1.3971 


1.9555 


1 


10 


272, 


48.7889 


3.5075 


.7339 


1 


10 


274 , 


49.3267 


3.3489 


.8322 


0 


9 


276, 


49,8658 


,4459 


.4249 


1 


10 


278, 


50.4053 


1.0606 


.6642 


1 


10 


280 . 


5q ,9456 


3.8362 


2.0458 


1 


10 


282, 


5l ,4871 


4.9476 


1.1952 


1 


10 


284, 


52.0292 


. 8078 


.5816 


1 


10 


286, 


52.5722 


1.3545 


. 63o6 


0 


9 


288, 


53.1158 


1.9762 


.9642 


1 


10 


290 . 


53.6600 


4.8418 


1 , 0683 


■ 0 


9 


292, 


54 , 2o52 


1.8616 


1.1776 


0 


8 


294 . 


54,7504 


2.3180 


1.2917 


0 


7 


296, 


55.2968 


1.1498 


.6087 


1 


9 


298 . 


55 . 844o 


2,4983 


.6028 


0 


8 


300 , 


56.3915 


1.6145 


.6383 


0 


7 
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FIGURE 7 
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NOISE- i 



RUN Jfl 
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FIGURE 8 
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FIGURE 9 
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FIGURE 10 
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FIGURE 11 
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FIGURE 13 




2Sa 



X-:>CfLl - 5.0CiE-fQJ UNJT3/JHCH. 

'f-3WLE - 100E-0Q 

•>*’•/ • 0 * I 'C » * ■* (J tl O .1 * '* » w V > <», /. -> ; 1* ^ , 5 

CENT2 EAST- SMOOTHING RESLl.TS, 



RUN lf\ 



tf 



*« 



( f f « Iff 

< « 

f f f f 

f « 

« f « f f f f 



i f 
4 4 i 



4 4 4 ( 
4 

4 4 f 
4 

4 14 4 



4 4 



« • 



^ • r 



(1 



FIGURE 14 
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where the rates of change of the direction cosines were high. As was 
expected the range smoothing results were the prime cause of the smoothed 
error . 

The system was able to smooth the tracks to within one and a half 
nautical miles inside 50 miles rather consistently, but it tended to respond 
rather strongly to the noise. It was noted, however, that this had the 
effect of causing no degradation in performance when the target began 
maneuvers. Also all the profiles yielded very similar smoothing results 
indicating that the smoothing constant values used were near optimum for 
the fixed constant, non -adaptive , system. 

In an attempt to more closely define the optimum values of the 
smoothing constants a program which varied the constants for each smoothed 
parameter from 0.05 to 1.00 in steps of 0.05 was written. This program 
is named ALPHA CHECK and listed in Appendix C . The run profiles used 
were similar to those listed earlier, but, the run time was shortened to 
100 seconds so that only the time that the target maneuvered was considered. 
The average mean square error the complete run was computed for each 
smoothed quantity over the range of the values of its smoothing constant. 

An autocorrelation analysis was also performed to evaluate this noise 
distribution'.' 'Tbocrosults are. shovm ?.n .Table IV.., Eor these runs the 
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AUTOCORRELATION ANALYSIS RESULTS OF PROGRAM 

ALPHA CHECK 
(51 SAMPLES) 





Mean 


Variance 


Largest Value of 
Autocorrelation 


RN 1 
Target 

Acfceleration 

Noise 


-0.1735 


0.6945 


0.3947 


RN 2 

Target 

Heading 

Acceleration 

Noise 


0.2439 


1.1572 


0.2662 


RN 3 

Target 

Pitch 

Acceleration 

Noise 


-0.0986 


.9622 


0.1774 


ATA 

Antenna 

Azimuth 

Noise 


0.264 


1.2022 


-0.2493 


EPSI 

Antenna 

Elevation 

Noise 


-0.0765 


0.8361 


0.2544 


R 

Radar 

Range 

Noise 


0.1438 


0.7524 


0.2783 


RDOT 

Radar 

Range 

Rate 

Noise 


-0.0661 


1.1695 


0.2805 


PD 

Radar 

Probability 
of Detectiqu > ; ’ • ■ 
Noise I I ’ 1 

r^i 1 — T^rrf 
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the probability of detection was one. This did not effect the results since 
no smoothing would have been performed when a contact was missed. The 
output would have been affected in an indeterminable way, however, by 
these predicted values. 

The results of these runs are shown in Table V. It was apparent 
from these results that there was considerable interplay between all the 
variables. Some rationale could be seen in this interaction, but it appeared 
nonlinear, since any change in the target's profile varied all the points of 
the optimum values for the smoothing constants . It was noted further that 
the knuckle in the curves at the optimum point was quite sharp indicating 
that system performance could be seriously affected when improper values 
of the smoothing constants were used. 

Due to the variations which were found in the optimum values of the 
smoothing constants for the chosen profiles it was not believed to be 
necessary to vary the parameters more widely since it was apparent from 
these runs that the system was very sensitive to changes in the profiles 
or the values of the system noise. 

The question must then be asked as to whether the re|}uirements 

i 

established for the analysis were too strict. In this regard specifically 
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OPTIMUM SMOOTHING CONSTANT VALUES 

FOR 

EACH TEST RUN 



Run 


Range 




Direction Cosines 


North 


East 


Down 


1 A 


0.1 


0„25 


0.25 


0.35 




0.05 


0„25 


0.25 


0.05 




0.15 


0.35 


0.35 


0.05 


2 A 


0o25 


0o90 


0.85 


0.65 


2 B 


LO 

o 

o 


0.65 


0.55 


0.10 


2 C 


0.35 


0.95 


0.95 


0.40 


3 A 


0.50 


0.60 


0.65 


0.55 


3 B 


0.45 


0 70 


0.55 


0.05 


3 C 


0.55 


0.75 


0.80 


0.05 
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the standard deviation of range noise and the 0.5 mile value chosen for 
the desired smoothing error become suspect since they were somewhat 
arbitrarily chosen. In rebuttal it is strongly felt, in the light of the high 
overall performance required of the entire system, its cost, and purpose, 
that the use of this relatively uncontrollable method of smoothing could be 
disastrous. Since the output of the smoothing section has considerable 
control over the performance of the system, specifically the accurate 
positioning of the missile, it is believed that a more adaptive system is 
desirable. It is further felt that a programmed logic system which would 
provide for varying the values of the constants would not be warranted 
since more sophisticated, and proven adaptive systems are available . 
Possible choices for such a method could be either the Least Squares , 
or Kalman filter methods . 

In short it must be concluded that the ^ Tracker 

system does not have the response or versatility required for a high 
performance air to air intercept problem. 
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APPENDIX A 

TABLE OF SYMBOLS USED IN COMPUTER PROGRAMS 



NOTE: If a symbol is used only in the ALPHA CHECK program it is so 

noted in parenthesis . 

ADEV: Absolute value of deviation obtained from RNDEV. 

This is used in the probability of detection computations. 



ADSMO: 

AESMO: 

ANSMO: 

AD: 

AE: 

AN: 

AR: 



Down, east, and north direction cosine differences 
from actual , and converted to radians . 



(ALPHA CHECK) Values of alphas for down, east, and 
north direction cosines, and for range, used for determining 
the optimum values . 



AERRM: Negative and positive desired limits on ADSMO, AESMO, 

AERRP: and ANSMO, a function of actual range. 
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ALPHAD: 

ALPHAE: 

ALPHAN: 

ALPHAR: 

ANOISE: 

AT: 

ATDOT: 

ATA: 
‘ATAE: , 
ATAD: 
ATADE: 

BETAD: 

BETAE: 

BETAN: 

BT: 

BTDOT: 



Position smoothing constants for down, east and north 
direction cosines, and for range. 



(ALPHA CHECK) Relative amount of radar parameter noise 
used in problem. 

Target heading and turn rate used in computation of target's 
position. 

Antennas actual and noisy azimuth angle measured positive 
clockwise , in radians and degrees . 



Rate smoothing constants for down, east, and north 
direction cosines. 



Target's pitch angle and pitch angle rate used in computation 
of targets position. 
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CIRERR- 


Three dimensional RMS actual and smoothed error difference 


CIRSMO: 


from actual target position. 


DEV: 


Random deviation obtained from random number generator. 


DT: 


Target position computation step size. 


DIRN: 


North direction cosine actual, noisy, smoothed and 


DIRNE: 

DIRNS: 

DIRNP: 


predicted . 


DIRE: 


East direction cosine actual, noisy, smoothed and 


DIREE: 

DIRES: 

DIREP: 


predicted. 


DIRD: 


Down direction cosine actual, noisy smoothed and 


DIRDE: 

DIRDS; 


predicted . 



DIRDP: 
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DIRND: 

DIRNDD: 


North direction cosine, rate and acceleration. 


DIRED: 

DIREDD: 


East direction cosine, rate and acceleration. 


DIRDD: 

DIRDD: 


Down direction cosine, rate and acceleration. 


DDSMO: 


(ALPHA CHECK) Down east and north direction cosine 


DESMO: 

DNSMO: 


smoothing differences from actual. 


DELTAT: 


Target decceleration time to final velocity at "1,75 g. 


DIVE: 


Target maximum dive angle, radians, corresponds to 60 deg 


EPSID: 


Antenna actual and noisy elevation angle measured positive 


EPSIDE: 

EPSI: 

EPSIE: 


downward, in radians and degrees. 





iimrrDSNTiAm 




» i»»> *i> ) >1) )i. >)> j 1 »»ti i)» ► » »> 

: 46 . of loo; : 



v) » ► tt 0 > ) » ) > t * * 



- r^* r. r r 



t > I 



y * * 



y y y 



y y y 



» » 1 * 



» (» (* (• (♦ 



ERAD: 


(ALPHA CHECK) Average mean square smoothed error for 


ERAE: 


down, east, and north direction cosines, and for range. 


ERAN; 


They are the quality indicators used for obtaining the 


ERAR: 


optimum values of the smoothing constants. 


F: 


A counter which fixes the target's pull out velocity. 


FPOINTSj 


(ALPHA CHECK) Number of radar returns received. 


GAMMAD; 


Acceleration smoothing constants for the down, east. 


GAMMAE: 


and north direction cosines. 



GAMMAN: 



GDN*; 


Target maneuver load factors for push over, pull out. 


GDP:* 


and turn. Note: GDN is a negative number. 



GT:* 



IT: 


Hollerith characters which identify the output plots . 


I: 


Number of frame times since problem start. 


K: 


Number of desired limits points plotted from problem start 



* Read into the program 
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NALPH: 


(ALPHA CHECK) Number of values of smoothing constants 
used in finding optimum value. 


ND: 


A zero or one depending on whether detection was missed 
or not. 


NERMAX: 


Number of points plotted indicating the desired limits on 
the smoothed errors. 


NGRAFl:* 


Setting NGRAFS to zero supresses various portions of the 


NGRAF2:* 


graph output. 



NGRAF3:* 



NOISE:* 


Number of times RNDEV is stepped prior to problem start to 
allow a different sequence of random numbers to be called. 


NPRINT: 


Number which separates pages of print out. 


NQ: 


Quality counter. 


NQC: 


Counter which reinitializes the initial conditions after the 
quality counter goes to zero. 



* Read into the program 
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NUMPTS: 



Number of curve points plotted . 



NUNIF: 

OM;* 

PD: 

PI: 

POVERD: 

POUTD: 

PSIO: 

PSIOD:* 

PSIT: 

PSITD:* 

PSITO: 

Q: 

* Read jp,to 



Counter in random number generator. 

Intercept 'r Mach number c 

Target detection probability for each frame time, a function 
of the actual range « 

3.14159265. 

Target pushover rate, rad/sec. 

Target pull out rate, rad/sec » 

Interceptor heading in radians and degrees . 



Target heading. 



Target initial heading, deg. and rad. Note a zero not an 



alphabetic o on target's initial positions. 



Number of frame. times since last radar return of range rate. 



the program . , , . , , , . ,,,,,,,, , , , , 
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R; 

RE: 

RNM: 

RNME: 


Actual and noisy range to target, ft and n.m. 


RERRM: 


Negative and positive desired limits on range, a 


RERRP: 


constant .5 n.m. 


RDOT: 


Actual and noisy range rate of target with respect to 


RDOTE: 


interceptor, positive opening, ft/sec. 


RDOTM: 


Range rate actual, noisy, predicted, and smoothed. 


RDOTME: 

RDOTMP; 

RDOTMS: 


n.m./seco 


RDDOTE: 


2 

Target noisy acceleration, n.m. /sec . 


RNl: 


Random accelerations of target for velocity, heading. 


RN2: 

RN3: 


and pitch. 
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RNMP: 

RNMS: 

RSMO: 

T: 

TA: 

TC: 

TD:* 

TM:* 

TS: 

TT I* 
VO; 

VT: 

VTDOT: 

VTF: 

* Read 



Predicted and smoothed range to target, nom. 



Smoothed range error from actual, n.m. 

Time since problem start. 

Total target iteration time since problem start. 

Iteration time which determines time for plotting 
desired limit points. 

Time to target dive . 

Target Mach number. 

Check time against iteration time to enter main program 
from target iteration . 

Time to target turn. 

Intercepter velocity, ft/sec. 

2 

Target velocity and acceleration, ft/sec. and ft /sec . 



Final target velocity, a function of the target's final 
altitude. 



into the program 
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XI: 

X2: 

XO: 

XONM: 

XTERR: 

XTO:* 

XT: 

XTE: 

XTONM: 

XTNME: 

XTNMS: 

XTSMO: 

Yl: 

Y2: 



Points which mark initial target location on the smoothed 
plot. 

Interceptor position along the north axis^ ft. and n.m. 



Difference between target's noisy and smoothed position, 

n.m. 



Target position along north axis, initial, actual, and noisy 
in feet, and initial, noisy, and smoothed position, n.m. 



Target's difference between actual and smoothed x position. 
Similar to XI and X2 . 
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YTERR; 


Similar to XTERR. 


YTO:* 

YT: 

YTE: 

YTONM: 

YTNME: 

YTNMS: 


Target’s position along the east axis, similar to X positions. 


YTSMO: 


Target’s difference between actual Y position and smoothed 



Y position. 



Zl: 

Z2: 


Similar to XI and X2 . 


ZETA: 


Target’s deflection angle from the horizontal, measured 



positive downward from the Intercepter^ rad. 



ZO:* 

ZONM; 


Interceptor's altitude, ft. and n.m. 


ZTERR: 


Similar to XTERR. 



Read into the program 
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ZTO;* 

ZT: 

ZTE: 

ZTNM: 

ZTNME: 

ZTNMS: 


Target’s altitude similar to X positions. 


ZTSMO: 


Target's difference between actual Z position and smoothed 



Z position. 



ZTF:* 


Final target altitude, ft. 


ZTPO: 


Target pull out altitude , ft . 



* Read into the program 
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APPENDIX B 

COMPUTER PROGRAM DETAILS 

A. Purpose: The program is designed to generate three dimensional 
target'dnterceptor geometries for use in analyzing optimum smoothing 
schemes . 

B. Limitations: The interceptor starts at the coordinate axis origin and 
may have any northerly heading. The target may be placed anywhere in 
the northern semi^circle with any initial heading. The positive axis 
directions are: X, north, Y, east, and Z, down, measured from the 
interceptc-r's altitude. The problem will stop anytime the target flies 
out of a 65 deg. cone from the interceptor's heading. 

C . Program Capabilities: The program will compute the actual range , 
range-rate, antenna azimuth and elevation for a radar frame time of 
two seconds. It will put noise on the above having the standard 
deviations of 2.5 n.m. for range, 8.72kts for range-'-rate , 0.745 deg. 

for azimuth, and 0.86 deg. for elevation. It also puts a haljf g acceleration 
noise on the target. It will compute the noitey target location, based on the 
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noise. It will further transfer the actual and noisy radar parameters 
into a moving north, east, down coordinate system which is fixed to 
move with the interceptor. It will determine whether detection was made 
or not each frame time, based on a detection probability curve, and 
perform a quality counter operation which prevents non-^detected targets 
from being smoothed after they are lost for a reasonable number of frames. 



D. Input data requirements: Fo|ur input cards are required. 

1. The first card reads in 16 Hollerith characters which identify the 
program on the print out and graphs . The format is 2A8 . 

2. The second card reads in the interceptor's altitude in feet, the 
target's X and Y coordinates in nautical miles, and altitude in feet, 
and interceptor and target heading in degrees. The format is 7F6.0. 

3. The third card reads in the target and the interceptor's Mach number, 
the target's maneuver g load capabilities for turn, pull out, and dive, 
(negative number) , and the problem time to target turn and/or dive. 

The format is 7F6.2. 

4 . The fourth card reads in a fixed point number that allows a different 
starting point for the subroutine RNDEV, and three fixed point ones or 
zeroes which control the graph output. The first, if zero, suppresses 
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the profile curves, the second, if zero, suppresses the smoothed 
difference plots„a.nd the third, if zefpo, suppresses the three dimensional 
rms smoothed difference plot. The format is 416, 



E . Requirements for reprogramming the smoothing section: The entering 
arguments of the smoothing and prediction section are RNME, DIRNE, 
DIREE, DIRDE, and NQC. The suggested statement numbers to use are 
the 600 series or above 2000. 



F. Quality indicators: The primary quality indicator is the three 
dimension rms smoothed error from the actual position. The three 
dimensional rms noise error is computed for comparison. The desired 
smoothing error on the individual quantities is 0.5 n.m. This could 
give a three dimensional rms smoothed error of as much as 0.875 n.m. 
The output curves have points plotted indicating these limits. 

G. Output: The print out contains the run identification read in, the 
problem time, the three dimensional -fms actual and smoothed error, 
whether detection was made or not, ND equals one or zero, and the 

- status of the quality counter, NQ. There are nine output curves. They 
plot the following: 
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1. The interceptor and target's actual and noisy positions in the 
north-east plane. 

2. The interceptor position and the target’s smoothed position in the 
north-east plane. 

3. The interceptor and target’s actual and noisy positions in the north- 
down plane. 

4. The interceptor’s position and the smoothed target’s position in the 
north-down plane. 

5. The range smoothed differences from actual and the desired limit 
points, versus time. 

6. The north direction cosine smoothed differences from actual and the 
desired limits points plotted as a function of actual range, versus time. 

7. The east direction cosine smoothed differences similar to the north 
direction cosine curve . 

8. The down direction cosine smoothed differences similar to the north 
direction cosine curve. 

9. The three dimensional f/ins smoothed error from actual and the 
desired limits points, vertsus time. 

H. Run time: The average run time is between 2 minutes and 30 seconds 

'*»»)'>> t ) t I I t . • » I ) I - » j * j ( » . » . ) ) ) ) 

and 2 minutes and! 5 Cv seconds . !’’! ! !, ! ! ! 

' ' > ‘ . j ) ) ) > ) ) ) I j > 

>> •>)>> • >1 J 
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PROGRAM THREE D 



THIS PROGRAM COMPUTES THREE DIMENSIONAL I NTERCEPTER- TARGET 
GEOMtiTRIER VOR THE"PUKPOSg' “Or; '01? T'rM'I;Z I ,VQ • V;AROUS PARAMETERS, 

THE RrOG<^^;M will t?.T.OP 'ANY’T.IME Th'AT ;THE'’TA;RGET FLIES OUTSIDE OF A 
C 65 DEGREE cb’N^ ' F'ROM ’ TWi^ rNTERt'EFTERS”HEAO TNG . 

C FOUR input cards ARE REQUIRED 

c 1 . The first card reads in i6 Hollerith characters which Ident- 

c IFY THE Program on the print out and graphs, the format is 2A8. 

C 2, THE SECOND CARD READS IN INTERCEPTERS ALTITUDE IN FEET, ' 

C targets X and Y COORDINATES IN NAUTICAL MILES AND Z COORDINATE 

C IN FEET, AND INTERCEPTER AND TARGET HEADING IN DEGREES. THE 

C format is 7F6.0. 

C 3. THE THIRD CARD READS IN TARGET INTERCEPTER MACH NUMBER, 

C target MANEUVER G LOAD CAPABILITIES FOR TURN PULL OUT AND DIVE 

C (NEGATIVE NUMBER), AND THE PROBLEM TIME TO TARGET TURN AND DIVE. 

C THE FORMAT IS 7F6.2. 

C 4. THE FOURTH CARD READS IN A FIXED POINT NUMBER THAT ALLOWS 

C A DIFFERENT STARTING POINT FOR THE SUBROUTINE RNDEV, AND THREE 

C FIXED POINT ONES OR ZEROS WHICH CONTROL THE GRAPH OUTPUT. THE 

C FIRST IF ZERO SUPPRESSES THE PROFILE CURVES, THE SECOND IF ZERO 

C SUPPRESSES The smoothed difference PLOTS, AND THE THIRD IF ZERO 

C SUPPRESSES The 3D RMS SMOOTHED ERROR PLOT, THE FORMAT IS 416. 



C DIMENSION AND READ STATEMENTS 



DIMENSION T(l51),X0NM(151) , YONMUSl) ,Z0NM(151),XTNM(15D,YTNM(151) 
1,ZTNM(151) ,RNM(151) ,ATAD(151),RNME(151) ,RD0TME(151) , ATADE(151) , 
2RNMS(l5l ) , RD0TM(l5l ) ,XTNME(l5l ) , YTNME(l5i ) , ZTNME(l5l ) , XTNMS ( l5i ) , 
3 YTNMS(i5i),ZTNMS( 151), I T ( 12 ) , D I RN ( l5i ) , D I RE ( l5l ) , D I RD ( l5l ) , 

4DIRNE(15i ) , DIREE dSl ),DIR0E(15 i ) , RDD 0TE(l5l ) ,RNMP(l5i ) ,RD0TMP(l5i ) , 

5, DIRNS(15i),DIRES(15i ),DIRDS(151),DIRND(151 ),DIRED(151),dIRDD(151) 

6, DIRNDD(i5i ) ,DIREDD(i5i) , dIRDDD(i5i) ,0IRNP(15i ) ,DIREP(i5i ) ,DIRDP 
7(i5i),EPSID(i5i),EPSIDE(i5i),RSM0(15i),ANSM0(i5i),AESM0(i5i),ADSM 0 
8(i5i), ZETAdSl) ,NQ(15 i ) ,XTSM0(i5i ) , YTSM0(l5l ) , ZTSMOdSl ) ,CIRERR(15 
9i) ,AERRP(3 q ), AERRM( 3o> ,RERRP(3 o) ,RERRM( 3o) ,TC(3 q) ,RERR(i5i) 

dimension XK2), Y1(2),Z1(2).X2(2),Y2(2),Z2(2),CERRP(30) 
l,XTERR(i5i ) , YTERRdSi ) ,ZTERR(l5i ) ,CIRSM0d5i ) ,N0d5j. ) - 

DO 2 J=i,6 

READ lOOdTd), 1=5,6) 

100 F0RMAT(2A8) 

READ ll,ZO,XT0NM, YToNM,ZTo,ZTF,PSIOD,PSITD,TM,OM,GT,GDP,GDN, 

ITT, TD, N0ISE,NGRAF1,NGRAF2.NQRAF3 

11 
C 

300 



FORMAT(7F6.0,/,7F6.2,/,4I6) 

INITIAL CONDITIONS 



Td) = 0 . 

TC(1)=0 . 

T A = 0 . !. >>'. 

TS = 0 . ’ ! ' 

NERMAX = 3o -- ■ 

PI=3.14i59265 
RAD=18o ./PI 
DI VE = -6o ./rad 
DT=.q 5 
NUMPTS=l5i 
NPRINT=68 



j * •» f » P 



NQC = 1 




NQ(1)=0 
NQ( 0 ) = 0 




XT0’=\f0NM*6b’76', 1 

YT0=YT0NM*6076,1 



XTsXTO 

YT=YTO 

ZT=ZTO 

PSlTo=PSlTD/RAO 
PSIO=PSIOD/RAD 
PSIT=PSITo 
AT=PSITO 
ATDOT = 0 , 

ATDDOTsQ . 

BT = 0 . 

8TDOT=0. 

BTDDOT = 0 . 

RN1 = 0 . 

RN2 = o , 

RN3 = o . 

NUNIF=122o7o3125 

C TARGET AND INTERCEPTER VELOCITY COMPUTATIONS 

VTDOT=0. 

IF(ZT0-37000. )63,69,69 
63 VT=TM*(1116.-.00398*ZT0) 

GO TO 57 
69 VT=TM*967,6 

57 IF(ZO-37000. >73,74,74 

73 VO = OM’»(1116.-. 00398*ZO) 

GO TO 75 

74 VO=OM*967.6 

75 IF(ZTF-10000 . >78,79,79 

78 VTF=1340 .+0 . 08*ZTF 
DELTAT=(VT-VTF)/56,4 
GO TO 58 

79 VTF=VT 
DELTAT=0. 

58 CONTINUE 

C this allows a VARIATION IN THEi RANDOM NUMBER SEQUENCE 

IF(NOISEU09,109,122 
122. DO 102 L=l, NOISE 
' 102 CALL RNDEV(NUNIF,DEV> 

109 CONTINUE 

C BEGIN MAIN PROGRAM DO LOOP 

DO 20 Is1%’15-2-'; ; •; , 



■ C 



CHECK’ I'VERATlOiV ’TINE ' ’A'GATNST ’RADAR fRAME TIME' 



22 IF (ABSF(TS-TA> -10 , £-5)25,25,24 



C 



CHECK MANEUVER TIMES AND/OR READ IN MANEUVER CONDITIONS 




62 of 100 



( r 

f f 



f r f 
f r 



c f 



. f 



« 











24 

39 

28 

29 

40 

41 

76 

77 

42 
46 

45 

34 

35 
203 

37 

200 

201 

43 



TA=TA*DT 

IF( TA-TT)39, 20^28 
ATDOT=0 . 

GO TO 29 

ATD0T=GT*32.i7/VT 

]ri^TA-TD.Mo*:4i,4i.^ 

l5TLJO.TiO.: ; 

Gb’ TO 42 ‘ * ' ” ’ • • 

IF(ZT-ZTF)76,77,77 
BTDOT=(GDP-l. )*32.17/VT 
GO TO 42 

BTDOT=(GDI^-l. )*32.17/VT 
VTDOT=-56.4 

IF( (PSITO+PI )-AT)46,46,45 
ATDOT=0 . 

AT=PSITo+PI 

IF(8T*DIVE)34,34,35 



> » I > > > 
) 

» > > 

» t » 
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BTDOT=0 . 

BT=DIVE 

IF (F) 2oi,2ol»2o3 

ZTP0=(1,-C0SF(BT) )*VT**2/( (GDP-1. )*32.17)+ZTF 
IF (ZTP0-ZT)38,37,37 

IF(F)201»201»200 

F=F-1. 

VTPO=VT 

BTDOT= (GDP-1. )*32,i7/VTP0 
IF(BT)38,43,43 

8TDOT=0 • 

BT = 0 . 



38 IF(ZT-ZTF)205.205,206 

205 ZT=ZTF 

206 IF(TA-(TD^DELTAT))51,51.2o2 
202 VTDOT=0. 

VT=VTF 



C ITERATION OF TARGETS POSITION 



5i VT = VT*VTD0T*DT-*-RN1*DT 

AT=AT+ATD0T*DT+RN2*DT**2/2, 

PSIT=AT 

BT=6T+BTD0T*DT+RN3*DT**2/2, 

C COMPUTE target POSITION 

XT = XT + VT*COSF(BT)*COSF( AT)*DT 
YT=YT+VT*C0SF(BT)*SINF( AT)*DT 
ZT = ZT-*-VT*SINF(BT)*DT 

C RETURN TO START TO CHECK AGAINST RADAR LOOK TIME 



GO TO 22 
25 TS = TS-^2, 



C 



59 



COMPUTE...INTERCE,PTERS. .POSITION, 



XO = VO*T(-I ) *OOSF-(,®'OIC'> 



AND A.nTU,A'. RAD, AR. .PARAMETERS 

<> ) I \> <i I 0 u I I 

0 J 0 I 0 ) 0 .) I 

1 I » I > > ) I I 

I > ) ) ) I > ) > I 



YO=VO*T( I )*SINF(PSI0) 

R = SQRTF( (XT-XO)**2-t-( YT-Y0)**2^(ZT-20)**2) 



EPSI=ASINF( (ZO-ZT)/R) 
IF(PSIO-PI )6o,59,59 
PS10=PSI0-2.*PI 



I 

I. « i 



• O C t I* t 



‘Am 



60 IF( YT-Y0)31»32,32 

31 ATA = -ACOSF( ( X T - XO ) / ( R*COSF ( EPS I ) ) ) -PS I 0 
GO TO 33 

32 ATA=AC0SF( (XT-X0)/(R*C0SF(EPS1)))-PSI0 




C 



" 'CHECK ’ irUARGei .rS’ .ST'I'L;:|*IN>'nADAR- 'B£A.H> IF NOT STOP 

‘ ■> O .1 . , , , . I > . ’ ’ 

’ , , • ’ ’ ’ ’ ’ ‘ • . I : ; ; ” 

33 1F( ATA-i’,i3‘^'5)^’6,26;5‘o' ’ ' •>' 



26 IF(ATA-^1. 1345)50,49,49 
49 IF(EPSI-i, 1345)48,48,50 
48 IF<EPSI+i;i345)50,27,27 
5q NUMPTS=(1-1) 

NERMAX=NUMPTS/10 
NPLOT=NUMPTS-NCTC 
IF(NUMPTS-60)15,i4,i 4 
15 NPRINT=NUMPTS 
l4 GO TO 1 

27 RDOT=-VT*COSF(PSIT-ACOSF( (XT-X0)/(R*C0SF(EPSI ) ) )*PI )*COSF(EPSI ) 
l-VO*COSF( ATA)*C0SF(EPSI ) 



C CONVERT ACTUAL RADAR PARAMETERS TO KTS, NM/SEC, DEG 



RNM( I )=R/6o76.1 
ATAD( I )5ATA*RAD 
EPSID( I )=EPSURAD 
RD0TM( I )=RDOT/6o76.1 
XONM( I )=XO/6o76.1 
YONM( I )=Y0/6q76.1 
ZONMC I )=Z0/6q76,1 
XTNM( I )=XT/6q76,1 
YTNM( I )=YT/6o 76.1 
ZTNMd ) = ZT/6o76.1 

C OBTAIN target MANEUVER NOISE 

CALL RNDEV(NUNIF,DEV) 

RN1=16.*DEV 

CALL RNDEV(NUNIF,DEV) 

RN2=i6 . /VT*DEV 
CALL RNDEV(NUNIF,DEV) 

RN3=i6,/vT*DEV 

C OBTAIN NOISY RADAR PARAMETERS AND CONVERT TO CONVENIENT UNITS 



C 



CALL RNDEV(NUNIF,DEV) 
ATAE=ATA+,o 13*DEV 
CALL RNDEV(NUNIF,DEV) 
RE=R+l5ooO .*OEV 
CALL RNDEV(NUNIF,DEV) 
EPSIE=EPSI+.o15*DEV 
CALL RNDEV(NUNIF,DEV) 
RD0TE = RD0T-^14.7i74*DEV 



ATADE( I )=ATAE*RAD 
EPSIDE( I )=EPSIE*RAD 



RNME( I ) = F}£/637.6ti. 
RD0TME( I|)=RD3 -t|E/60 7'6* 




't i> 
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COMPUTE NOISY TARGET SPACE STABILIZED POSITION 



IF(PSIO-l, )7o,7o,7l 

70 XTE=XOfRE*COSF( ATAE+PSI0)*C0SF(EPSIE) 
YTE = Y0^RE*SINF(ATAE-'-PSI0)*C0SF(EPSIE) 






( f < 



c < 



c 

c 



71 

72 



GO TO 72 

XT£=XO+Rb*COSF( ATAE-PSlO)*COSF(EPSIE) 
YTE = YO-^RE*SInF(ATAE-PSIO)*COSF(EPSIE) 
2TE=Z0“RE*SINF(EPSIE) 




XTNME( J )=XTE/6076.1 
Y'sNKEU >=,V7£/;6o76;l 
Z'^NKEl ri= 2 TE /^60 74.J: 

* * * • > 0 » ♦ j ^ 'k > 
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COMPUTE ACTUAL DIRECTION COSINES FOR INTERCEPTER WITH ROLL 
AND PITCH ZERO 



DIRN{ I ) = COSF(ATA )*C0SF(PSI0)-SINF (ATA)*SINF(PSIO) 

DIRE( I )=COSF( ATA)*SINF(PSI0)-*-SINF(ATA)*C0SF(PSI0) 

DIRD( I ) = -SrNF(EPSI ) 

C COMPUTE NOISY DIRECTION COSINES 

DIRNE( I )=C0SF(ATA£)*C0SF(PSI0)-SINF(ATAE)*SINF(PSI0) 

DIREE( I )=C0SF(ATAE)*SINF(PSI0)-*-SINF( ATAE)*C0SF(PSI0) 

DIRDEI I )=-SINF(EPSIE) 

C OBTAIN DETECTION PROBABILITY AND DETERMINE IF TARGET WAS 

C DETECTED 



CALL RNDEV (NUNIF.DEV) 

IF(RNM( I )-2o. >21,23,23 
21 PD = 2,36-, 043*RNM( I) 

GOTO 52 

23 IF (RNM( I )-5o . >36,47,47 

36 PD=1,95«.0239*RNM( I ) 

GO TO 52 

47 IF (RNM( I >-8q . >65,66,66 

65 PD=1.55-.016*RNM( I > 

GO TO 52 

66 IF (RNM( I >-110 , >62,60,60 

62 PD=.88-,0077*RNM(I> 

GO TO 52 

68 PD= ,25-.00i9*RNM( I > 

52 CONTINUE 
ADEV=ABSF(DEV> 

IF(PD-ADEV>53»54,54 

C quality COUNTER COMPUTATIONS 

53 ND( I >=0 

NQ( I >=NQ( I-1>-1 

IF (NQU > )25o, 151,152 ,, 

250 NQ(I>=0 
GO TO 151 

54 ND( I )=1 

NO( I >=NQ( I-l>+2 
IF (10-NQ{ I > >25i,i53,l53 

251 NQ( I >=10 
GO TO 153 

c shoo>t>hin:g aVid PREiiiciTioVi se-:i:io’n 

* > >> > )»>))> 
f)IT) llrtO !>•) kOVJ i>) » >)» ) ) ) I 

C no detection conditions 

C IF QUALITY COUNTER IS ZERO PRINT THE FOLLOWING AND START NEXT 

C FRAME TIME 



I 



4 



151 CIRERR( I )=999.9999 



'WSBSESS^SSSSIB^ 



CIRSMO( I )=999.9999 




XTNMS( I ) =XTNM( I ) 
YTNMS( I ) ?YTNM( I ) 
ZTNMS( I )5Z7NM( I ) 
ND( I ) = 99 ,* /*’ , 



C 

c 



IF QUALITY COUNTER IS NON ZERO ERROR VALUES EQUAL PREDICTED 
ONES, COMPUTE ONLY NEXT PREDICTED VALUES 



152 RNME( n=RNMP( I ) 

RDOTMEC I )=RDOTMP( I ) ' 

DIRNE( I )=DIRNP( I ) 

DIREE( I )=DIREP( I ) 

DIRDE( I )=DIRDP( I ) 

GO TO 61 

C DETECTION CONDITIONS 

153 if (N0C)6i,6i,i54 

C initial contact conditions at start or AFTER CONTACT HAS 

C BEEN LOST AND REGAINED 

154 NOCsNOC-l 
RNMP( I )=RNME ( I ) 

RDOTMPC I )=RDOTME( I ) 

DIRNP( I )=DIRNE( I ) 

DIREP( I )=DIREE( I ) 

DIRDP( I )=DIRDE( I ) 

RDDOTE( I-1)=0 . 

DIRND( I-1)=0 . 

DIRNDD( I-1)=0. 

DIREDC I-1)=0 . 

DIREDDC I-1)=0 . 

DIHDD( I-1)=0 . 

DIRDDD(I-1)=0. 

61 CONTINUE 

C SMOOTHING CONSTANT COMPUTATIONS 

ALPHAR=.45 

ALPHAN=.43 

8ETAN=ALPHAN**2/(2 . -ALPHAN) 
GAMMAN=2.*BErAN/ALPHAN-(ALPHAN+8ETAN ) 

ALPHAE=ALPHAN 
ALPHAD= .35 

8ETAD=ALPHAD**2/(2 , -ALPHAD) 
GAMMAD=2.*BETAD/ALPHAD-(ALPHAD+8ETAD ) 

C SMOOTHING AND PREDICTION COMPUTATIONS 



C 



RANGE 



■) ) ) o 



RNMS( I )=RNMPt I*H-'ALPHARi(FiNM5^ I I.3-RNMP(I ) ) ’ ’ ‘I , 

RDDOTEC I ) =RDDOTE( I ’1 ) +1 . / ( 2 , *Q ) ★ ( RDOTME ( I ) -RDOTMP ( I ) ) 
RNMP( I+1)=RNMS( I )*2.*RD0TME( I ) •'■2 . *RDDOTE ( I ) 

RDOTMP( l+D=RDOTME( I ) +2 , *RDDOTE ( I ) 



C 



NORTH DIRECTION COSINE 




• X' 



« I • 
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D1RNS( I )=DIRNP( I )-t-ALPHAN*(DlRNE( I)-DIRNP( I ) ) 

DIRND( I )=DIRND( I -1 ) -^2 . *D I RNDD ( I -1 ) ••■BE T AN/2 . * ( D I RNE ( 1 ) -D 1 RNP ( I ) ) 
DIRNDD( I )=DIRNDD( I-1)+2.*GAMMAN/4.*(D1RNE( I )-DIRNP( I ) ) 

DIRNP( 1*1 ) = DIRNS( I )+2 .*DIRND( I )*2.*DIRNDD( I ) 

C . ;'EA3:r DIRECTION CC3IN2 i ’ : ’ : 

' ' I'O , , , / • 

DIRES( I )=DIREP( I )*ALPHAN*(DIREE( I )-DIREP( I ) ), 

DIRED( I )=DIRED{ I-1)*2.*DIREDD( I-1)*BETAN/2,*(DIREE(I )-DIREP< I )) 
DIREDD(I)=DIREDD( I-1)*2.*GAMMAN/4.*(DIREE( I )-DIR£P{ I ) ) 

DIREP( I*1)=DIRES( I )*2.*DIRED( I ) *2 . *D I REDD ( I ) ! 

C DOWN DIRECTION COSINE 

DIRDS( I )=DIRDP( I )*ALPHAD*(DIRDE( I )-DIRDP( I ) ) 

DIRDD( I )=DIRDD( I-1)*2.*DIRDDD( I-1)*BETAD/2.*(DIRDE( I )-DlRDP( I ) ) 
DIRDDD( I )=DIRDDD( I-1)*2.*GAMMAD/4,*(DIRDE( I )-DIRDP( I ) ) 

DIRDP( 1*1 )sDIRDS( I )*2 ,*DIRDD( I )*2 .*DIRDDD( I ) 

C COMPUTE SMOOTHED TARGET POSITION 

ZETA( I )=AC0SF(DIRDS( I ) )-PI/2, 

XTNMS( I )=XONM( I )*RNMS( I ) *COSF ( ZET A ( I ) )*DIRNS( I ) 

YTNMS( I )=YONM( I )*RNMS( I ) *COSF ( ZET A ( I ) )*DIRES( I ) 

ZTNMS( I ) = ZONM( I )*DIRDS( I )*RNMS( I ) 

C COMPUTE smoothed AND ERROR DIFFERENCES FROM ACTUAL 



XTSMO( I )=XTNM( I )-XTNMS( I ) 

YTSMO( I )=YTNM( I )-YTNMS( I ) 

ZTSMO( I )=ZTNM( I )-ZTNMS( I ) 

XTERR( I )=XTNM( I )-XTNME( I ) 

YTERR( I )=YTNM( I )-YTNME( I ) 

ZTERR( I )=ZTNM( I )-ZTNME( I ) 

RSMO( I )sRNM( I )-RNMS( I ) 

C COMPUTE THE ACTUAL AND SMOOTHED 3d RMS ERROR VALUES 



CIRERRI I )=SQRTF(XTERR( I )**2+YTERR( I )**2*ZTERR( I )** 2 ) 

CIRSMOI I )=SQRTF(XTSMO( I ) **2*YTSM0 ( I ) **2-^ZTSMO ( I ) **2 ) 

C CHECK TO SEE IF THE ACOSF ARGUMENT HAS BEEN EXCEEDED# 

C IF IT HAS SET THE COSINE EQUAL TO ONE AND PRINT OUT THE VALUE 



IF (DIRNS( I )-l. )2ll,2ll,2l0 



210 PRINT 212, (T( I ) ,DIRNS( I ) ) 

212 F0RMAT(6X2HTs F8 . 4 , 2 X6HD I RNS = F8 , 4 ) 

DIRNS( I )=1. 0 

211 IF(DIRES( I )-i, )22l, 221,220 



220 PRINT222, (T( I ),DIRES( I ) ) 

222 F0RMAT(6X2HT=F8.4,2X6HDIRESsF 8.4) 
DIRES( I )=l,o 



221 

230 

232 



I F ( D I RDS I ) -1, 0 ) ?3l , 231,, 2?0 
PRINT23^, (;T('I-5,DIRD£U 5) 



} ) > 
* •) 



> > 
} ) ) > 



FORMAT ( 6X2H7sffF8; 4, 2X6KC'i RD3sF3 . 4 ) - 



DIRDS( I )=1,0 
231 CONTINUE 



• 



C 



CONVERT THE DIRECTION COSINES TO DEGREES. FOR GRAPH OUTPUT 



I • * 
« % 



• • • 
• • 



ANSM0(I)=AC0SF(DIRN(I))-AC0SF(DIRNS(I)) 
AESMO( I ) = ACOSF(DIRE( I ) ) - ACOSF ( D I RES ( I ) ) 
ADSMO( I )=ACOSF(DIRD( I ) ) - ACOSF ( D I RDS ( I ) ) 




C 



COMPUTE THE DESIRED LIMIT POINTS FOR THE GRAPH OUTPUTS 



IF.<Tf;(K/ fT( I ) U20»l-’0»l2l 

120 TC(K‘-*-i")='rcvvo vix’. •' * ■’ -- 

AERRP(K)=ATANF( ,5/RNM( I ) ) 

AERRM<K)'s-AERRP(K) 

RERRP(K)? ,5 
RERRM(K)=-.5 
CEHRP(K)=.875 
K = K*1 

121 CONTINUE 

C IF DETECTION WAS MISSED INCREASE Q IN THE ROOT SMOOTHING 

C COMPUTATIONS 

IF (ND( I ) >55,55,56 

55 Q = Q-*-l. 

GO TO 2o 

56 Q=l, 

2o T( I+1)=T( I)’*-2, 

C END OF MAIN PROGRAM DO LOOP 

1 CONTINUE 

C OUTPUT CONTROL 

C COMPUTE markers WHidH INDICATE THE TARGETS INITIAL POSITION 

C ON THE SMOOTHED PROFILE CURVES - 

X1(1) = XT0NM-*-.5 

XI (2 ) = XT0NM- .5 

Y1 (1) = YT0NM-*-.5 

Y1(2 )sYT0NM-.5 

Zl(l)=ZTO/6o76,i+.o5 

Z1(2)=2To/6o76,i-.o5 

X2(1)=XT0NM+.5 

X2(2)=XToNM-.5 

Y2(1)=YToNM-.5 

Y2(2)=YToNM+.5 

Z2(1)=ZT0/6076,i- ,q5 

Z2(2)=2To/6o76.i*,o5 

C GRAPH STATEMENTS '' 

IF (NGRAF1)27o,27o,27i 
- 271 CONTINUE 

DO 101 I=l»4 
101 IT(I)=8H 

IT(D = 8HGENTZ NE 
IT(2)=8H, plane a 

IT(3)=8HCTUA,L.i* . 



IT(4)=8HNCISE 
1 = 1 



CALL DRAWINUMPTS, YONM,XONM,1,0#^H , I T , 15 , , i5 , , q , 3 , 2 , 2 , 6 , 6 , 1 
ILAST) 

CALL. DRAW(NUMPTS,YTNM,XTNM,2,0,4H ) 

CALL DRAW(NUMPTS, YTNME,XTNME,3,o,4h ) 




4 4 
4 4 4 



• • • 

4 

• • 



1 



DO 103 1=1,4 

103 ITU)=8H 
IT(1)=8H0ENTZ NE 
IT(2)=8H PLANE S 
IT(3)=8HM00THED 

IT(4.)=8HP0?,I,T,10N 

I =1 r ' r ' ' 

CAl I.'- r.RAW'vNliMPTS,'MDNM’>’xONMil,C,^H 
ILAST ) 

CALL DRAW(NUMPTS, YTNMS,XTNMS,2,0,4H 
CALL DRAl^ (2, Y1,X1,2,0,4H ) 

CALL DRAW ( 2 , Y2 , X2 , 3 , Q » ) 

DO 104 1=1,4 

104 IT(I)=8H 
IT(1)=8HGENTZ ND 
IT(2)=8H PLANE A 
IT(3)=8HCTUAL ♦ 

IT(4)=8HN0ISE 

1=1 

CALL DRAW(NUMPTS,XONM,ZONM,1,0 
ILAST) 

CALL DRAW(NUMPTS,XTNM,ZTNM, 2,0,4H 
CALL DRAW(NUMPTS, XTNME , ZTNME , 3 , 0 » 4H 
DO 110 1=1,4 
110 IT(I)=8H 

IT(1)=8HGENTZ ND 
IT(2)=8H plane S 



» ®»,l5, ,0,3f2,2,6f6i^l, 

i 



,IT,15,,1,.5,o,0,2,2^6^6,1 , 

) 

) 



IT(3)=8HM00THED 

IT(4)=8HP0SITI0N 



CALL DRAW(NUMPTS,XONM,ZONM,1,0»4H , I T , i5 . , 1 . 5 , q , 0 » 2 , 2, 6, 6 , 1 , 
ILAST ) 

CALL DRAW(NUMPTS,XTNMS,ZTNHS,2,0,4H ' ) 

CALL DRAw(2,X1,Z1,2,o,4H ) 

CALL DRAW ( 2, X2, Z2. 3, 0 , 4H ) 

27o CONTINUE 

IF (NGRAF2)281,281,282 
282 CONTINUE 

DO 1q5 1=1,4 

105 IT(I)=8H 
IT(1)=8HGENT2 RA 
IT(2)=8HNG6 SMOO 
IT(3)=8HTHING RE 
I T(4 )=8HSULTS 
1=1 

CALL DRAW(NUMPTS,T,RSMO,1,0,4H , I T , 5q , , 1 . , 3 , Q , 2 , 2 , 6 , 6 , 1 , LAST ) 

CALL DRAW(NERMAX,TC,RERRP,2,2,4H ) 

CALL DRAW(NERMAX,TC,RERRM,3,2,4H ) 

DO 106 1=1,4 

106 IT(I)=8H 



IT(1)=8HGENTZ NO 
IT(2)=8HRTH SMOO 
IT(3)=8HTHING RE 
IT(4)=8HSULTS 

1=1, ■ ‘ ' 

CALL L'RAV! 1 NUMPTS> T j Al^SMO, 1 , 0^ Wirl 

CALL’ DRAW vNE'rK'A'X', 7S’,A2r'r?',’2, Z,'4H 
CALL DRAW(NERMAX,TC,AERRM,3,2,4h 



, II', "^0 .’*'»0*L,3^o,2,2,6,6j1,LAST) 
) 



DO 116 1=1,4 
116 IT(I)=8H 



IT(1)=8HGENTZ EA 
IT(2)=8HST Smoot 




« 1 



IT(3)=8HHING RES 
IT(4)=8HULTS 
1 = 1 

.CALL DRAW(NUMPTS.T,A£SMO,1,0,4H 
CALL- ;O.RAW ( NcPMAX, f C # AERRp , 2 / 2, 
CALL DRAW'{NEPMAx:TC,A'ERRrt,3,2,“H’ 
do’ 108 1 = 1,4 ' ' ’ 

108 IT(I)=8H 

IT(1)=8HGENTZ DO 
IT(2)=8HWN SMOOT 
IT(3)=8H^ING res 
IT(4)=8HULTS 
1=1 

CALL DRAW(NUMPTS,T,ADSMO,1,0,4H 
CALL DRAW(NERMAX,TC,AERRP,2,2,4H 
CALL DRAW<NERMAX,TC,AERRM,3,2,^H 

281 CONTINUE 

IF (NGRAF3)283,283,284 

284 CONTINUE 

DO 285 1=1,4 

285 IT(I)=8H 
IT(1)=8HGENTZ 3D 
IT(2)=8H RMS SMO 
IT(3)=8H0THING R 
IT(4)=8HESULTS 

1 = 1 

CALL DRAW(NUMPTS,T,CIRSMO,1,0,'»H 
ILAST) 

CALL DRAW(NERMAX,TC,CERRP,3,2,4H 
283 CONTINUE 

C PRINT STATEMENTS 




, lT,5o ,01»3,0,2,2,6,6,1,LAST) 

i ) 

' ) ' : : 



r 



, lT,5p,i «01>3,0«2,2,6,64l,LAST) 

) 

) 



I 



,IT,5o*, <5, 0 , 0 , 2, 2, 6, 6 ^ 1 , 

) 



I 

1 



PRINT 95 

95 F0RMAT(1H1,27X17HSM00THING RESULTS) 

PRINT 64,(IT(I),I=5,6) 

64 F0RMAT(/,32X,2A8) , 

PRINT 107 

107 FORMaT( /,6X66HT = TIME, R NM = ACTUAL RANGE, C I RERR = 3D I MENS I ONAL RMS 
1 ACTUAL ERROR, , / , 6X64HC I RSM0=3 DIMENSIONAL RMS SMOOTHED ERROR i 

2, ND = 1 DETECTION, ND = 0 N0,/,6X 68HDE TECT I ON , NQ = QUALITY COUNTER, W 
3HEN NO=o target CONSIDERED LOST AND ,/,6x24HALL NINES ARE PRINTED. 

4 , I 

PRINT 96 I 

96 F0RMAT( /,i5xiHT,4x4hr NM, 7x6hC I RERR , 4x6hC I RSMO, 2X2HND, , 

12X2HNQ) 

PRINT 97, (T( I ) ,RNM( I ) ,CIRERR( I ) »CIRSMO( I ) ,NDU )#NQ( I) , I=1,NUMPTS) 

97 format ( 7X,1F10.0#3F10,4,2I4) 

2 CONTINUE 

400 STOP 
END 
END 





i o > 



t t 



t 0 






# > .' , ♦ ) . > ) ) 

» r nr * ) #.) > >» i 

'" ’’ ' ' ’ 'Apt>endix’ C ■ '> 

Computer Program of Program ALPHA CHECK 



Only the smoothing and prediction and output sections are shown 
since the preceeding portions are the same as in THREE D. 



A V 
I » 



' . .71 cf iroo: 

\ . . n . 









i * 



47 IF (RNM( I )-80. )65,66,66 

65 PD = 1,55-. 016*RNM( I ) 

GO TO 52 

66 IF (RNM( I )-110 , )62,68,68 
62 PDs. 88-, 0077*RNM( I ) 

GO .TO 52 , 

68 PD= . 254,J)o19*RNH( 1 )^. 

52 COMTINOjE 
AD6V=ABSF(DEV) 

ADEVsO , ‘ 

IF(PD-ADEV)53»54,54 

53 NO( I ) = 0 
GO TO 20 

54 NQ(I)=1 

20 T( I+1)=T( I )^2. 

1 CONTINUE 




) > # i 






I 



C SMOOTHING COMPUTATIONS 



FPTSsNUMPTS 

NALPHS20 

c range smoothing constant analysis 



ALPHARs , 05 
DO 301 J=1 ,NaLPH 
P = l. 

AR( J)=ALPHAR 

ERAR( J)*o • 

DO 3q5 I=1,NUMPTS 
IF (NQ( I ) )3o2,3o2,3q3 
302 IF(P)304,304,305 
3q4 RNME( I)sRNMP( I) 

RDOTMEI I )sRDOTMP( I ) 

GO TO 3o6 

3q3 IF(P)3o6,3o6,3o8 

308 P=P-1. 

CTC = T( I )/2.^l, 

RNMP( I ) sRNME( I ) 

RD0TMP( I )=RD0TME( I ) 

RDD0TE( I-1)=0 . 

3o 6 RDD0TE( I )=RDD0TE( I-1)+1./(2.*Q)*(RD0TME( I )-RD0TMP( I ) ) - 
RNMS( I)sRNMP( I )+ALPHAR*(RNME( I )-RNMP( I ) ) 

RNMP( I*1)sRNMS( I ) +2.* ROOT ME ( I )+2,*RDD0T£( I ) 

RDOTMP( I-^l)*RDOTM£( I )f2.*RDD0TE( I ) 

RSMO( I )=(RNM( I )-RNMS( I ) )**2 
ERAR( J)sERAR( J)*RSMO( I ) / ( FPTS-CTC ) 

IF(NQ(D) 55,55,56 

55 Q=Q^1. 

GO TO 3q5 

56 Q=l, 

3o5 CONTINUE 



301 ALPHAR = ALPHAR-*-, 05 
1=1 



P R ^ N T > 11 4 ) > » ) > i > > 

ll4 FCRMrAT'!i:HlLl7y,o6'H:0PT,l’MUMi SMOo’thVnG CONSTANT ANAL YS I S , // , 6X62hOUAL 

ITY {NDiCATC'Ri TliE- MINI MUM AVERAGE MEAN SQUARE ERROR ) 



PRINT 604, ( ( IT( I) » 1=5,6) , ANOISE) 

604 FORMAT(//,20X,2A8,2X7HANOISE=F4.2) 
PRINT 115 

ll5 F0RMAT(//,29 Xi6HRANGE, ALPHAR ) 

PRINT 116, (AR(J),ERAR( J),J=1,NALPH) 







I • « ( 



* « 



* * 



^rr 



116 FORMAK//, ( (23X7HALPHAR= F A . 2 , 2X5HERAR= F6.4)/)) 




C NORTH DIRECTION COSINE SMOOTHING CONSTANT ANALYSIS 



ALPHAN=.o5 

DO 32i„ lt = l,NALPH. . ' , 

a I >«■' i* • 

AN(K)-ALPHAN ' ' 

ERAN(K)=0 . 

BETANsALPHAN*w2/(2.-ALPHAN) 

GAMMANS2 . *BET AN/ALPHAN- ( ALPHAN+BETAN ) 

DO 325 I=l,NUMPTS 
IF(NQ( I ) 1322,322,323 

322 IF(P)324,32A,325 

324 DIRNE( I ) =DIRNP( I ) 

GO TO 326 

323 IF(P)326,326,328 
328 P=P-1. 

DIRNP( I )=DIRNE( I ) 

DIRND( I-1)=0 . 

DIRNDD( I -l)sO . 

326 DIRND( I )=DIRND( I-1)-*-2.*DIRNDD( I -1 ) +BETAN/2 . * ( D I RNE ( I )-DlRNP( I ) ) 
DIRNDD( I )=DIRNDD( I -1 ) +2 . *GAMMAN/4 . * ( D I RNE ( I )-DIRNP< I ) ) 

DIRNS( I ) =DIRNP( I )+ALPHAN*(DIRNE( I )-DlRiNP( I ) ) 

DIRNP( I+1) = DIRNS( I )+2.*DlRND( I ) +2 . *D I RNDD ( I ) 

DNSMO( I )=(DIRN< I )-DIRNS( I ) )**2 
ERAN(K)=ERAN(K)*DNSM0( I )/(FPTS-CTC) 

325 CONTINUE 

321 ALPHAN=ALPHAN+, 05 
1=1 

PRINT 117 

117 FORMAT(//,2ox31HNORTH DIRECTION COSINE, ALPHAN ) 

PRINT 118, ( AN(K),£RAN(K) ,K=1,NALPH) - 

118 FORMAT(//, ( (23X7HALPHAN= F4 . 2 , 2X5HERAN= F10.8)/)) 



C EAST DIRECTION COSINE SMOOTHING CONSTANT ANALYSIS 



ALPHaE= , 05 
DO 33l L = 3.,NALPH 
P = l. 

AE(L) =ALPHAE 

BETAE=ALPHAE**2/(2.-ALPHAE) 

GAMMAE = 2.*BETAE/ALPHAE-(ALPHAE-^BETAE) 
ERAE(L)=0 , 

DO 335 I=1,NUMPTS 
IF (NQ( I ) >332,332,333 

332 IF(P)334,334,335 
334 DIREP( I)=DIREE( I) 

GO TO 336 

333 IF(P)336,336,338 



- 338 P=P-1. 

DIREP( I)=DIREE( I) 

DIRED(I-1)=0. 

DIREDD( I-l)so . 

336 DIREDA I /=JDIRED( ;-1//+2 . -d:rEDD( I '1 ;+SCTAE/2. i(DlREE( I )-OlREP( I ) ) 
DIRED'J{ r)'=DlREDt-( I-l)’-^2.*6Af=^MAE/4 *(DtPEE( :0-DIREP( I ) > 

DIRES( i)^DlREP{ l)*MLrHAE*(DIREE( I ;-DiREP( I / ) 

DIREP( I+1)=DIRES( I )-*-2.*DIRED( I ) *2 . *D I REDD ( I > 

DESMO( I)=(DIR£( I )-DIRES( I) )**2 
ERAE(L)==ERAE(L)+DESM0( I )/(FPTS-CTC) 

335 continue 



331 ALPHAE=ALPHAE*, 05 




APPENDIX C 



1=1 

PRINT 119 

119 FORMATdHl, 20X3 ih EAST DIRECTION COSINE, 
print 161, (AE(L),ERAE(L),L=1,NALPH) 

161 FORMAT!//, ( (23X7HALPHAE= F4 . 2 , 2X5HERAE= 




AUPHAE ) 

FIO .8)/) ) 



I •» » % > » 



, DOWi^>,.D.i RECT I 0\^ COiJlNE' £'MOO]TH;I'NG ^OfiSTA^T .'ANALYSIS 

' j» o a 

ALPHAD=.05 



» » > •» rt • 



» » > \ > 1 



DO 341 M=1,NALPH 
P = l. 

AD(M)=ALPHAD 

BEl AD = ALPHADh.*2/(2, -alPHAD) 
GAMMAD=2.*8ETAD/ALPHAD-(ALPHAD+8ETAD ) 

ERAD(M)=o, 

DO 345 I=l,NUMPTS 
IF(NQ( I ) >342,342,343 
342 IF(P)344,344,345 
344 DIRDE( I )=DIRDP( I ) 

GO TO 346 



343 IF(P)346,346,348 
348 P=P-i, 

DIRDP( I)=DIRDE( I ) 
DIRDD( I-1)=0 . 



DIRDDD! I-1)=0 , 

34 6 DIRDD( I )=DIRDD( I-1)-*-2.*DIRDDD( I-1)-*-8ETAD/2.*(DIRDE(I )-DlRDP( I ) ) 
DIRDDD(I)=DIRDDD(I-1)^2.*GAMMAD/4.*{DIRDE(I)-DIRDP(I)) 
DIRDS(I)=DIRDP(I )-*-ALPHAD*(DIRDE( I)-DIRDP(I)) 

DIRDP( I'^1) = DIRDS( I )+2.*dIRDD( I ) ■*•2 . *0 I RDDD ( I > 

DDSMO( I )=(DIRD( I )-DIRDS( I ) )**2 
ERAD(M)=ERAD(M)+DDSM0(I)/(FPTS-CTC) 

345 CONTINUE 

341 ALPHAD = ALPHAD>^. 05 



1=1 

. PRINT 121 

121 FORMAT!//, 20X31HDOWN DIRECTION COSINE, ALPHAD ) 
PRINT 123, ! AD!M) ,ERAD!M) ,M=1,NALPH) 

123 FORMAT!//, ! !23X7HALPHAD= F4 . 2, 2X5HERAD= FlO.8)/)) 
IF !NGRAPH)27 o, 270,271 
271 CONTINUE 

DO llo 1=1,4 
110 IT!I)=8H 

IT!1)=8HGENTZ AL 
IT!2)=8HPHAR VS 
IT!3)=8HAVG mean 



IT!4)=8HSQ ERROR 
1=1 

CALL ORAW!NALPH,AR,ERAR,o,0»4H 
DO 111 1=1,4 
111 IT!I)=8H 



,IT,,2,o,,0,0«2,2,5,6,i,LAST) 



IT!1)=8HGENTZ AL 
IT!2)=8HPHAN VS 
IT!3)=8HAVG mean 
IT!4)=8HSQ error 

1=1 

CALL DR'\W(NALPH,'ANtS.RAN,0>0,4P’ 
DO 112 i =1, 4 ■ 

112 IT!I)=8H 

IT!1)=8HGENTZ AL 
IT!2)=8HPHAE vs 
IT!3)=8HAVG mean 
IT!4)=8HSQ error 



: ’ > I T , .2, c . >*0 >0*# 24 2, 5 , 6 , 1 , LAST ) 

' ' >, 



< f 



( c < 



t t 



t t t t t f 



f f * 



APPENDIX C 



113 



1=1 

CALL DRAW(NAlPH,AE,ERAE»0»0»^H » I T , . 2 , q . * 0 * 0 
DO 113 1=1, A 

I T •( I ) ='8 H* ••>• •« >«>»•• I 

iTiiii=:8nqEiNTz: AL : :•* : !’* i : 

I T (*2 ) i-8HPHAl? vs '•* ’ 



2,2»5»6»1»LAST) 



IT(3)=8HAVG mean 
IT(4)=8HS0 error 
1=1 

CALL DRAW(NALPH,AD,ERAO,0,0,^H , I T , . 2 , 0 . , Q , 0 * 2 , 2 , S , 6 , 1 , L^ST ) 

270 CONTINUE ' 

AN0ISE=AN0ISE^.5 
2 CONTINUE 
END 
END 



r Ok : 



i « « 

i 

• * \ 



< V* • 



t ^ * 



» » » 



App'ejndiJ?.* C 

Results of Runs 2A and 3A 



. 7^, of. i.on . « M. j . 



i • 



( { 



% « « 



± 1jJ-i V X 



SMOOTHING RESULTS 
RUN 2A 




T=TIME. R NM= 


ACTUAL RANGE, CIRERR 


=3DIMENSI0NAL 


RMS ACTUAL ERROR, 


CIRSMO-3 .D.IMENSiONAL R^.S SMOOTHES 


. €RR-ORvND-«l .DETECTION, ND»0 NO 


DETECT'! OM,: N3 


i;QUA4I TY; 


.C.GUNTpR; *.WH 


EN;N-Qso: TA41GE;T 


pONSIDERED LOST AND 


ALL NINE’S 'ARE 


PRi NTE-D*. 


• »'•» • 




f 




T 


R NM 


CIRERR 


CIRSMO 


NO 


NQ 


. I 

• 


14 . 15q8 


.2429 


.2429 


1 


2 


2. 'i 


14,5361 


1.8578 


.1884 


0 


1 




14,9288 


5.1681 


2.3338 


1 


3 


6. . 


15,3280 


2.6577 


2.4972 


1 


5 


8. 


15.7357 


.2391 


1.3032 


1 


7 


10 . 


16 ,l53o 


3.2430 


2.1581 


1 


9 


12. 


16,5815 


.5384 


1.3833 


1 


10 


lA. 


17.0281 


2.8164 


.5864 


1 


10 


16. 


17,4867 


.2840 


.4138 


1 


10 


18. 


17.9446 


3.1760 


1.2653 


1 


10 


20. 


18 , 4i55 


.3837 


.6776 


1 


10 


22. 


18,8975 


1.3631 


1 .0167 


1 


10 


24. 


19.3867 


.26o7 


.4953 


1 


10 


26. 


19.8780 


.3462 


.2626 


1 


10 


28. 


2o .3618 


1.0^72 


.4439 


1 


10 


30. 


20.8512 


.4633 


s .4537 


1 


10 


32. 


21.3431 


.8665 


.65q7 


1 


10 


34. 


21.8341 


3.2o78 


1.2478 


1 


10 


36. 


22.3337 


2.3855 


.5537 


1 


10 


38. 


22.8417 


.6766 


.2386 


1 


10 


40 . 


23.3525 


2.2829 


1.o736 


1 


10 


42. 


23.8669 


.4162 


.7837 


1 


10 


44, 


24,3896 


.8367 


, 49o5 


1 


10 


46. 


24.9098 


2.7580 


1.1853 


1 


10 


48. 


25,4187 


1,8427 


.4945 


1 


10 


50 . 


25.92q7 


1.8184 


.8741 


1 


10 


52. 


26,4281 


2.0666 


1.4061 


1 


10 


54. 


26,9446 


.1707 


.7198 


1 


10 


56. 


27.4663 


.3516 


.3665 


1 


10 


58. 


27,9841 


2.4726 


1.2136 


1 


10 


60 . 


28,4961 


2.8985 


1.9486 


1 


10 


62. 


29.0163 


.9102 


.8421 


1 


10 


64. 


29.54o5 


2.4250 


1.5345 


1 


10 


66. 


30 .0547 


.2888 


.7757 


1 


10 


68. 


30.5636 


2.8572 


1.0178 


1 


10 


70. 


31.0732 


1.1682 


. 7910 


1 


10. 


72. 


31.5752 


2.5874 


.9229 


1 


10 


74. 


32.0719 


4.1386 


.9257 


0 


9 


76. 


32.5737 


4.3305 


.9190 


0 


8 


78. 


33.0744 


.7393 


.2930 


1 


10 


80 . 


33.5744 


.4837 


.2606 


1 


10 


82. 


34.0735 


3.1434 


1.3395 


1 


10 


84. 


34.5622 


6.2619 


2.0776 


1 


10 


86.",: 


•55'. G’347 


; '■X'.37'08/' 


:3..,64 99; 


'■•'3, : 


U 


88’. j 


55 -.4 993 


1,39*2^5 •' 


; .15853; 


i; : 


•1'3 : : 


90.” ■ 




■^’.^2'91’” 


' i“.l751' 


^”1 ' 


10 ‘ ’ 


92. 


36.4471 


.8413 


.9765 


1 


10 


94. 


36 . 92o5 


4.1630 


1.0l70 


0 


9 


96, 


37.4q29 


.7908 


1.q545 


0 


8 


98, 


37.8898 


2.9077 


1.8114 


1 


10 


100. 


38.3598 


2.2822 


1.9747 


1 


10 - 
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102 . 


38.8208 


104, 


39.2771 


106. 


39 . 7219 


108. 


40.1699 


110 . 


40.6274 


U2, . 


.41. ,0.9 7 5 


1-14; ; 


’.:4;i.58Sl 


i’i.6’. : 




118, 


42.5788 


120 . 


43.0949 


122. 


43,6149 


124. 


44.1265 


126, 


44.6289 


128. 


45.1289 


130 , 


45.6308 


132. 


46.1254 


134. 


46.6149 


136, 


47.1069 


138. 


47.6015 


140 . 


48.1068 


142. 


48.6124 


144, 


49.1099 


146. 


49.6123 


148. 


50.1305 


l5o . 


5q .6560 


152. 


51.1413 


154. 


51.5500 


156. 


51.8594 


158. 


52.0482 


160 . 


52.1140 


162. 


52.0591 


164. 


51.9006 


166. 


51.6697 


168 . 


51.3834 


170 . 


5l .0395 


172. 


50.6556 


174. 


50,2446 


176. 


49.8208 


178. 


49,3989 


180 . 


48.9821 


182. 


48.5692 


184, 


48.1119 


186. 


47.5735 


188, 


46.9689 


190 . 


46.3166 


192. 


45.6382 


194. 


44.9558 


196. 


44.2740 


198. 


43.5920 


200 . 


42,9103 


202. 


42.2292 


204. 


41.5480 


206. 


40.8671 


208. 


40.1867 


2lo . 


39. 5,0 63 


212. 


: 3:8.A'?5 9; 


214. 


’ 2;8.-l45f^ 


2l6. 


37.4664 


2l8. 


36.7869 


220 . 


36.1085 


222. 


35.43o3 


224, 


34.7529 



.2352 


1.6414 


1 


10 


. 0081 


1.6821 


0 


9 


.3226 


1.7492 


0 


8 


.5017 


.6377 


1 


10 


.5599 
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